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Abstract: The move towards simplified models for Run II of the LHC will allow for 
stronger and more robust constraints on the dark sector. However there already exists 
a wealth of Run I data which should not be ignored in the run-up to Run II. Here we 
reinterpret public constraints on generic beyond-standard-model cross sections to place 
new constraints on a simplified model. We make use of an ATLAS search in the mono jet + 
missing energy channel to constrain a representative simplified model with the dark matter 
coupling to an axial-vector Z'. We scan the entire parameter space of our chosen model to 
set the strongest current collider constraints on our model using the full 20.3 fb -1 ATLAS 8 
TeV dataset and provide predictions for constraints that can be set with 20 fb -1 of 14 TeV 
data. Our technique can also be used for the interpretation of Run II data and provides a 
broad benchmark for comparing future constraints on simplified models. 
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1 Introduction 

In recent years Effective Field Theories (EFTs) have become a popular framework with 
which to constrain the dark sector at the LHC [1-12]. In the simplest cases, the dark cou¬ 
plings and mediator masses are combined into a single effective energy scale, A, 1 leaving 
this and the dark matter mass, moM, as the only free parameters for each effective oper¬ 
ator. EFT constraints have the advantage of being relatively model-independent, allowing 
constraints to be placed across a broad range of models and parameters. In addition they 
facilitate an easy comparison with direct detection experiments via the shared energy scale 
A. However it is now clear that EFTs must be used with extreme care at LHC energies, 
where the energy scale is large enough that the approximations used in the construction of 
EFTs can not be assumed to be valid. At these energies and luminosities, the energy carried 
by the mediator is usually larger than the mediator mass, violating the EFT approxima¬ 
tions, except in the case of large mediator masses or for dark-sector couplings approaching 
the perturbativity limit [11-21]. Depending on the mass and width of the mediator, this 
can lead to EFT constraints that are either stronger or weaker than the constraints would 
be on a UV-complete model, reducing their utility and making their validity questionable. 

1 Also called M* in the literature. 
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One solution is to rescale EFT constraints, by truncating the simulated signal such 
that only events for which the EFT approximation are valid are used to derive constraints 
[14, 22, 23]. This weakens constraints but at the same time makes them substantially more 
robust , which is critical when considering bounds on beyond-standard-model parameters. 
Whilst this technique has the advantage of maintaining some of the elegance of EFTs, 
it also has the serious disadvantage that it does not make full use of all potential signal 
events available in a UV complete model and so does not address the region of parameter 
space where EFT constraints are too weak. To constrain this region we need to consider 
models where the mediator can be resolved. On the other hand, the parameter space of 
full, well-motivated models such as supersymmetry [24] or extra dimensions [25] is broad, 
and by focusing solely on such models we run the risk of missing more generic signatures 
of the dark sector. 

Hence, the usage of simplified models is now advocated by a number of groups [26-31]. 
Here we will use publicly available ATLAS constraints on the monojet + missing energy 
channel to constrain a simplified model with dark matter coupling to the standard model 
via exchange of an axial-vector Z' mediator. The original search was used to constrain 
EFTs, however the same data and analysis can be used to constrain a simplified model of 
choice through the model-independent limit on the visible cross section contribution from 
beyond-standard-model processes. Such a reanalysis only requires simulation of the signal 
in the new model for each point in parameter space. 

Simplified models have the advantage of a relatively small set of free parameters, and 
do not encounter the same validity problems as EFTs. However, the parameter space is still 
larger than for EFTs, which often necessitates arbitrary choices for one or more parameters 
in order to constrain the remaining free parameters. Here we will instead leave the dark 
matter mass, mediator mass, and coupling strength all as free parameters which we scan 
over and constrain in contours. We derive limits using the full ATLAS 8 TeV dataset [23] 
and predict the limits that can be set with 20 fb -1 of 14 TeV data. 

In Section 2, we outline the choice of simplified model that we will be constraining. In 
Section 3, we describe our technique for converting the model-independent constraints on 
the visible monojet cross section into constraints on this simplified model. In Sections 4 
and 5 we present our results, before we give our concluding remarks in Section 6. The 
Appendix contains validations of our limit-setting and various approximations used, and a 
discussion of the relic density calculation. 

2 Model 

We consider a widely-used benchmark simplified model where Dirac DM interacts with 
the SM via a Z'~ type mediator. This is described by the following Lagrangian interaction 
term: 

£ = - E - S^75)<?] - Z'p [xt m 0dm - 0dmT5)x], (2- 1 ) 

/ 

where g(, gf are respectively the vector and axial-vector coupling strengths between the 
mediator and quarks (i = q) and DM (i =DM). This is a well-motivated simplified model 
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that has been studied extensively, including searches by CMS [32] and ATLAS [33], and 
numerous other groups both in the UV complete and EFT limits, e.g. Ref. [28, 30, 34]. It 
is part of the wider family of dark Z* portal models which have been studied previously 
in e.g. [35-39]. The LHC is relatively insensitive to the mixture of Vector/Axial-vector 
couplings [23], however this ratio has a large effect on the sensitivity of direct detection 
experiments to this model. A vector coupling induces a spin-independent (SI) WIMP- 
nucleon scattering rate, while an axial-vector coupling induces a spin-dependent (SD) rate 
[40]. Current bounds on SI interactions are much stronger than those on SD, to the point 
where direct detection constraints are generally stronger than LHC constraints on models 
with pure vector couplings, and vice-versa for pure axial-vector couplings, as seen in e.g. 
Ref. [41]. For this reason we consider a pure axial-vector coupling, setting g^ M — g^ — 0, 
and defining $dm = #dM’ 9q — 9^- ^ or s i m phcity, as is common, we assume that the quark- 
mediator coupling g q is the same for each species of quark. We require that ^dm 5 gq < 47r 
individually in order for the couplings to remain in the perturbative regime, however in 
practice limiting Tos/M ends up constraining our parameter space further. 

For the model we consider, the total on-shell width of the axial-vector mediator is 
given by: 


Tos 


9 dm M ( 1 ^ m DM ) 3 ^ 2 


12tt 


©( Af — 2771 dm) T 


E 

q 


g q M( 1 — 477i^/M 2 ) 3 / 2 

47r 


0(M — 2771 




( 2 . 2 ) 

where M is the mediator mass. With the assumption that g q is equal for each flavor of quark 
this width can become very large, for example rising above Tos M at g q = g D M ~ 1-45 
when g q = ^dm- This width assumes no additional decay channels aside from quarks and 
DM, however it is conceivable that the mediator could decay to standard model leptons or 
additional dark sector particles. Given that the the structure of a possible dark sector is 
unknown and that couplings to leptons are more appropriately constrained by searches for 
dilepton resonances [42, 43], we confine ourselves to the more ‘minimal’ model where the 
mediator couples only to quarks and DM. For a study of how the limits change when the 
width is manually made larger (without considering specific additional decay modes) see 

[30]. 

In the event generation we will use a Breit-Wigner form for the Z' propagator with 
constant on-shell widths: 

1 


A z ,(s,M,T) = 


s — AL 2 + i ALTos 


(2.3) 


This is questionable from a first-principle Quantum Field Theory perspective [44, 45]: 
it amounts to a replacement of sT(s)/M T(s)M in the imaginary part of the Dyson 
resummed self-energy, which clearly can only be motivated for small widths for which 
s ~ M 2 . Additionally we identify T(s) as its on-shell value Tos; This is formally correct 
in our tree-level context, but can only be expected to be a good approximation when 
Tos/AT <C 1. A completely consistent theoretical treatment of our model in the part of 
parameter space where we can’t assume that Tos/M <C 1 hence requires a more careful 
consideration of the propagator than we employ in our simulation, and we therefore use 
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cross section reweighting for the parts of parameter space where we are not conservative 
to correct for this. In Figure 1 we show a comparison of line shapes for various values 
of Tqs/M with fixed M, mo m using the naive Breit-Wigner shape we employ and a fully 
kinetic propagator, which reflects that our simulation will become increasingly poor as 
Tos/M increases, as shown previously in [35]. There are also further corrections which 
can not be included by using a kinetic width at tree level as we do here, but these are 
subleading and this treatment should be sufficient for our purposes. 

We have estimated the ratio of final cross sections using a Breit-Wigner propagator 
and a kinetic one by convolving them with parton density functions (PDFs) and requiring 
that we can produce a XX pair on-shell and found that within this approximation the 
choice of propagator makes for modest differences for most of our parameter space, with 
large effects in specific regions: for ttidm M and Tqs/M > 0.5 our propagator gives a 
cross section that can be as much as 50% too low compared to a kinetic propagator. In 
the other extreme end where ttidm > M and Tos/M > 0.5 we see an opposite effect where 
our cross section can be several times too large. Since our propagator in this region is not 
conservative we will reweight the visible cross section we get at the end of our simulation 
whenever ^Breit-Wigner > & kinetic using the ratio of the two (where these cross sections are 
estimated as described above) to allow our limits to remain robust. 



Figure 1: Line shapes for our Breit-Wigner propagator 2.3 (blue) and a kinetic propagator with the 
replacement TosM sT(s)/M in the imaginary part of the self-energy (red) for various values of Tos/M. 
The normalisation is arbitrary and differs between the plots to allow for a straightforward shape comparison, 
and both axes are linear. M — 500 GeV and ttidm — 100 GeV which corresponds to a typical parameter 
space point. 
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3 Reinterpreting Monojet Constraints 

Our signal prediction is obtained by implementing the model in the FeynRules [46] and 
MadGraph5_aMC@NLO 2.1.2 [47] framework to generate leading order (LO) parton 
level events using the NNPDF2.3 LO PDFs [48]. These are matched to Pythia 8.185 [49] 
using the MLM algorithm with a matching scale of 80 GeV 2 for showering and hadronisation 
using tune 4C. We generate XX + 0, 1, and 2 jets in the matrix element before matching 
to the parton shower. We use the default MadEvent factorization and renormalization 
scales which in this case both are approximately the transverse mass of the XX 

system. Our approach only makes leading order + parton shower (LOPS) predictions 
compared to the next-to-leading order + parton shower (NLOPS) predictions used in a 
similar study of CMS results [50] in [28], which means we suffer from larger theoretical 
uncertainty due to scale dependencies which we can attempt to estimate by varying our 
choice of jur : f by a factor of two. This shows a weak dependence on the choice of scales of 
IJJj* for a few representative choices of M, ttidm which is clearly not a realistic estimate of 
the uncertainty: previous studies [51-53] with other choices of scales have found fixed-order 
NLO corrections ranging from ^ 20 — 40%. We do however note that based on the results 
in [53], we expect fixed-order NLO corrections to ultimately be modest after matching to 
a parton shower and applying the ATLAS monojet analysis cuts since the parton shower 
dilutes differences, helped by the loose cuts on additional jets. As such they should have a 
limited impact on our quantitative results and be negligible for qualitative results. 

We analyze the generated events using the ATOM framework [54, 55] based on Rivet 
[56]. We first divide the final state into topological clusters and find jets with the anti -kt 
algorithm [57] using R = 0.4 in Fast Jet [58]. We then perform a smearing of the px of 
these jets based on typical values for the ATLAS detector, leaving the E^f lss unsmeared 3 . 
Our procedure has been validated by recreating the limits set on an EFT operator by 
ATLAS and the results of this validation is given in Appendix A. 

Some past constraints on simplified models have used a fixed benchmark width. In 
this case, the cross section is only sensitive to the product $dm * 9 q and not to the couplings 
individually; Further, this easily factorises out, 

d<r{gDM,g q ) = (#dm • g q ) 2 d(j(gBM = g q = 1), (3.1) 

which simplifies the analysis since the coupling affects only the magnitude of the signal, 
not the spectral shape. Including the physical width complicates things, since now both 
the magnitude and signal spectrum have a dependence on both $dm • g q and g q /gr>M- This 
results in necessary complications if one wants to present 2D contour limits on $dm • g q 
when the width is known. However it is possible to make an approximation for the cross 
section in the resonant region as a oc p 2 p^ M /ros (for fixed M, ttidm) and a oc g^g^ M for 
the off-shell region, which allows us to set limits on $dm * g q while avoiding a scan in this 

2 Chosen to correspond to the matching scale used in the original ATLAS EFT interpretation. 

3 We are not aware of any ATLAS E™ 1SS smearing values which could be unambiguously applied to our 
case, based on the results in [59] we expect the plateau to have been reached for all our signal regions 
however. 
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dimension, leaving only M — tudm as free parameters for any given choice of g q /gr>M- This 
approximation should work well for the part of parameter space where Tos M [17, 60], 
and we present a full study of this approximation (including the effect of the mediator 
shape reweighting) in Appendix B which further motivates restricting the parameter space 
to Tos/M < 0.5. 

We also include relic density constraints by finding out which parts of our parameter 
space would result in a larger relic abundance than observed experimentally. The details 
of this calculation are described in Appendix C. 

4 20.3 fb -1 8 TeV Limits 

We make use of 90% C.L. limits set using the full ATLAS 8 TeV dataset [23] which should 
allow us to set the strongest collider constraints yet on our model due to improvements in 
the selection and the use of more signal regions compared to the full dataset analysis by 
CMS [61]. The analysis defines signal regions based on E^f lss and initially only requires the 
leading jet > 120 GeV with the additional requirement that 2 -p^ > E™ 1SS for the harder 
signal regions, and does not veto events due to additional jets as long as A</>(jet, E™ 1SS ) > 
1.0. The signal regions are given in Table 1. 


Signal Region 

SRI 

SR2 

SR3 

SR4 

SR5 

SR6 

SR7 

SR8 

SR9 

#y ss & 2 • $ > [GeV] 

150 

200 

250 

300 

350 

400 

500 

600 

700 

ATLAS cr^f 0 CL [fb] 

599 

158 

74 

38 

17 

10 

6.0 

3.2 

2.9 


Table 1: Signal region definitions in the full dataset 8 TeV analysis and 90 % CL exclusion limits on the 
visible cross section from BSM contributions. 


4.1 Discussion of constraints 

The results of our full scan are shown in Figure 2. They cover more parameter space and are 
typically stronger than previous constraints using CMS results, although our limits tend to 
be slightly weaker in the off-shell region which can be traced back to our very conservative 
mediator shape reweighting. We also make an attempt at properly mapping out the limits 
on the g q /gr>M < 1 part of parameter space for the first time. We see that we can exclude 
mediator masses up to ~ 1.3 TeV for g q — (/dm = 1 and ^ 0.9 TeV for g q — (/dm = 0.5 
at 90% C.L. Lower g q /gr*M values are relatively well-constrained by monojet searches as 
expected due to resonant enhancement of the cross section. Additionally because g q is 
small, dijet constraints are relatively weak for this part of parameter space as we will 
discuss below. Although briefly studied in [28, 35], to our knowledge these results are the 
first to map out the limits in this part of parameter space in detail with the full dependence 
on M and tudm- 

The area inside the black line in Figure 2 indicates the region where the constraint 
on ^/g q gDM is stronger than the coupling strength which gives the correct relic abundance. 
Naively, M — tudm values inside the region outlined by the black line would therefore lead 
to a larger relic abundance than observed. However, if the DM is not produced thermally 
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Figure 2: Current ATLAS monojet constraints at 90% C.L. on our model with the full 8 TeV dataset. 
The dashed white line shows where Tos/M = 0.05. The black dots are interpolation knots in M — ttidm 
space. The region inside the black line is naively ruled out by relic density contraints. For details see text. 


or if the DM has other annihilation channels not considered here (e.g. to leptons), this 
constraint is relaxed. Therefore this contour does not strictly rule out any region, but 
rather indicates that the ‘simplest’ models of DM are expected to lie outside this contour. 
Limitations aside, this is a useful way to compare collider measurements to cosmological 
observations, and can be elegantly implemented in a simplified model context as we show 
here. 

4.2 Limits from dijet resonance searches 

Since our Z' couples to quarks one could attempt to make use of limits from dijet resonance 
searches [32, 62-64] to further constrain the model as done in for example [35]. This 
certainly gives much stronger constraints for the g q /gr>M > 1 part of our parameter space 
which is why we don’t study this part of parameter space (we show results for g q /gr>M = 2 
in our validation study in Appendix B.l). For lower ratios of couplings we need to be 
careful since dijet resonance searches generally assume narrow resonances and in this part 
of parameter space the dark sector branching starts contributing to the width considerably. 
The dashed white line on the plots show where the width of the mediator becomes narrow 
enough to potentially violate such constraints assuming no additional dark sector decays 
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(we take this to be Tqs/M < 0.05 to be conservative, but note that there are recent 
searches [64] which have constrained much wider resonances). Comparing to the detailed 
Z' dijet analysis in [65] and the recent ATLAS update in [62] and assuming the results 
won’t change drastically when using an axial-vector coupling compared to a vector one, 
we see (note the difference of a factor of 6 in the definition of g q ) that the values of g q for 
which Tqs/M < 0.05 for g q /gr> m = 1/2,1/5 in our model typically are smaller than the 
values currently constrained by dijet searches, but g q /gr>M — 1 might be better constrained 
by dijet searches in the part of parameter space inside the dashed white line. Realistically 
the dijet resonance constraints will apply for much wider resonances than we allow here 
once you perform a proper analysis instead of relying on constraints set assuming narrow 
widths, as done in [35] at parton level which suggests that monojet searches give the 
strongest constraints for light mediators but dijet searches are more constraining for heavy 
mediators even for low values of g q / gr>M- Due to the sensitive dependence on the width it is 
however worth stressing that since the model we assume here has no additional dark sector 
or standard model decays for the mediator, constraints set on dark matter mediators with 
dijet resonance searches in this part of parameter space can not be considered conservative: 
the width we use is the minimum width assuming equal coupling to each generation of 
quarks, and small changes to the dark sector can make a large difference to this width. 
This problem is not as pronounced when g q /gr>M > 1 since the width then is dominated by 
SM decays, which further motivates using dijet constraints over mono jet ones in this part 
of parameter space. We also note that since the width can be large, interference effects 
with Z/y* should be properly taken into account when using dijet searches to constrain 
these models - we expect interference to play a similar role as in Drell-Yan [66] and have 
checked that this appears to be the case but a detailed analysis is outside the scope of this 
paper. 

It is also possible to make use of dijet angular distributions which are sensitive to 
wider resonances than the dijet mass spectrum [67, 68] and therefore can be considered 
more robust than dijet resonance contraints. As shown in the parton level study in [35] 
these can also be competetive with resonance searches in some parts of parameter space, 
but we won’t consider them further here. 

5 14 TeV Predictions 

We make use of the public results in [22] as estimations for the expected backgrounds 
and hence expected cross section limits at 14 TeV with 20 fb -1 of data from the ATLAS 
detector assuming a low average number of pile-up collisions (/i = 60) 4 . The details of 
the simulation and analysis are largely the same so we assume the detector performance 
will not degrade which is well-motivated under the current upgrade predictions, but our 
analysis is changed to mirror that in [22]: we use a constant leading jet pt cut of 300 GeV 
and the signal regions have been redefined as detailed in Table 2. The estimated reach 


4 We thank David Salek for providing us with the exact numbers for the background estimates and 
expected limits. 



with 20 fb 1 of 14 TeV data for g q /gr>M = 1/5,1/2,1 is presented in Figure 3. The black 
line again indicates the correct relic density as discussed in Section 4.1. 


Signal Region 

SRI 

SR2 

SR3 

E™ ss > [GeV] 
ATLAS <j“ p - 95% CL [fb] 

400 

28 

600 

4.5 

800 

1.5 


Table 2: Signal region definitions in the 14 TeV analysis and expected 95 % CL exclusion limits on the 
visible cross section from BSM contributions. 




Figure 3: Predictions for the reach of the ATLAS experiment at 95% C.L. with 20fb _1 of 14 TeV data 
with ii — 60. The dashed white line shows where Tos/M = 0.05. The black dots are interpolation knots 
in M — ttidm space. The region inside the black line could naively be ruled out by relic density contraints. 
For details see text. 


6 Conclusion 

As the LHC approaches Run II there is a clear move towards supplementing EFT analyses 
with simplified models, as a stronger and more robust way to constrain the dark sector. 
These same arguments apply to Run I data, and thus it is useful to reinterpret existing 
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constraints on the dark sector in the simplified model framework. This has the added ben¬ 
efit of allowing clearer benchmarks and comparisons for future studies of simplified models 
at higher LHC energies and luminosities. We have demonstrated this with constraints 
on a simple Z ' model, with an axial-vector coupling. The parameter space for simplified 
models spans a minimum of 4 dimensions, making the parameter scan and visualisation of 
the subsequent constraints more challenging than for EFTs. The common restriction to 
2-D slices of parameter space does allow for easy comparison between several constraints, 
but reduces our knowledge of the model as a whole. Here we instead scan over the full 
4-D parameter space, presenting results as contours, allowing us to retain the maximum 
information possible on constraints of a dark sector. By making some well-motivated ap¬ 
proximations we can perform such scans in an accurate and computationally reasonable 
way. 

This allows for a more complete understanding of the strengths of the monojet channel 
for constraining dark sectors and facilitates comparison to contraints from other experi¬ 
ments and astronomical observations, as shown here by comparing to dijet and relic density 
constraints. Whilst the scope of this analysis is limited to a single simplified model, this 
technique shows good prospects for the (re)interpretation of constraints across a broader 
model-space. 
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A Validation of Procedure using EFT Limits 

We use the ATLAS monojet analysis for 10.5 fb -1 of 8 TeV data [69] to validate our 
procedure and approximations, so we apply the following cuts: We require at most two 
jets with px > 30 GeV and \rj\ < 4.5, with |t 7 j1 | < 2 and A0(j2, E™ 1SS ) > 0.5 where j\ and 
j 2 are the leading and subleading jet respectively. We define four signal regions based on 
pt p and E™ 1SS , given in Table 3. 

Our overall limit-setting procedure has been validated by recreating the ATLAS limits 
set on A for the D8 EFT operator which corresponds to our simplified model. A comparison 
for SR3 is presented in Table 4. We consistently overestimate the limit by a few percent 
which reflects the less advanced nature of our detector simulation, however the agreement 
is good enough for our purposes as we have sub-2% differences for ttidm values which are 
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Signal Region 

SRI 

SR2 

SR3 

SR4 

pf & £y ss > [GeV] 

120 

220 

350 

500 

ATLAS al jf° CL [fb] 

2800 

160 

50 

20 


Table 3: Signal region definitions in the 10.5 fb 1 8 TeV analysis and ATLAS 95 % CL exclusion limits 
on the visible cross section from BSM contributions. 


relevant for us. Note that we only perform the comparison for SR3 as it usually is the most 
discerning signal region and the only one for which ATLAS results are reported, however 
we assume the results are similar for the other signal regions. Similarly we assume this 
agreement carries over to our analyses of the full 8 TeV dataset and our 14 TeV predictions, 
which is well motivated since the full dataset 8 TeV analysis was conducted under similar 
conditions and due to the stated ATLAS upgrade goals for the upcoming higher energy 
LHC run respectively. 


m D M [GeV] 

ATLAS 95% CL on A [GeV] 

Our 95% CL on A [GeV] 

Difference [%] 

<80 

687 

700 

+1.9 

400 

515 

525 

+1.9 

1000 

240 

250 

+4.2 


Table 4: Comparison of limits set on the D8 EFT operator by ATLAS [69] and us using only SR3. 


B Validation of Cross Section Reweighting 

Our limits set using results from Ref. [69] using interpolation in M — moM — #dm * 9q are 
presented in Figure 4, limits set using the cross section approximation including the width 
mentioned in Section B.l are presented in Figure 5, and the ratios of the limits set in the 
two cases are presented in Figure 6. To visualise the breakdown of our approximations we 
extend the limit of our parameter space to Tos/M < 1. 

Values of g q /gr>M > 1 are hardly probed at all by monojet searches as evident from our 
results for g q /gvM = 2: such models are much better constrained by dijet searches which 
motives not including these in our main study. 

B.l Using a cross section approximation including the width 

We compare our results to ones obtained by reweighting the cross section for a single value 
of <7dm • g q to see how well the simple cross section approximations: 

a oc { »«4 M/r ° S if M > 2 "‘ DM (B.l) 

f 9q9DM li M < 2m DM 

reproduce the full results. Additionally we perform a separate reweighting to correct for the 
Breit-Wigner shape of the propagator as for the full results, although only before finding 
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Figure 4: Our results using interpolation in M — tudm — 5dm • g q space. The dashed white line shows 
where Tos /M = 0.05. The black dots are interpolation knots in M — ttidm space. See the text for further 
details. 


the limit using B.l 5 . The results using this reweighting are presented in figure 5 and the 
ratio of the limits in Figure 6. As expected they work well for sufficiently small values 
of 9dm • 9q that Tos M. As Tos/M -T 1 the cross section approximations become 
inaccurate and the less accurate shape reweighting becomes an issue. This motivates us to 
use these cross section approximations for our study of the full 8 TeV ATLAS dataset and 
our 14 TeV predictions while remaining robust by further limiting the parameter space to 
Tqs /M < 0.5. 


C Comparison with Relic Density Constraints 

If dark matter was produced thermally in the early universe, there is a simple relationship 
between the thermally averaged dark matter self-annihilation cross section (crp) ann , and 
the observed relic abundance Udm^ 2 - For a given model, this allows us to find the cou¬ 
pling strength which provides the correct relic abundance as a function of (m x ,M). This 

5 This could be improved by using an iterative approach but we avoid doing so here to keep the procedure 
simple. 
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Figure 5: Our results using interpolation in M — tudm space and the cross section approximations detailed 
in the text. The dashed white line shows where Tos/M = 0.05. The black dots are interpolation knots in 
M — ttidm space. See the text for further details. 


scenario is by no means a certainty; If the observed dark matter was produced through 
some mechanism other than thermal production, or if some new physics has an effect on 
the connection between the self-annihilation rate and the abundance at freezeout, this rela¬ 
tionship breaks down. At the same time, thermal dark matter is a well-motivated scenario, 
and is a useful way to get a sense of the regions of parameter space within which we expect 
the gravitationally-observed DM to lie. 

The observed relic abundance can be approximated as [70] 


Ddm^ 2 — 


2 x 2.4 x 10- 10 GeV -2 

( ( 7f’)ann 


(C.l) 


Combined with Planck constraints of D^/i 2 = 0.1199 ± 0.0027 [71], we see that (crp) ann ~ 
4.0 x 10 -9 GeV -2 for thermal relic DM. Rather than using this approximation, we use 
the well-known formalism described in, for example, Refs. [72, 73] to constrain (crp) ann by 
simultaneously solving an expression for the freezeout temperature as a function of (crp) ann , 
and the relic abundance as a function of both (crp) ann and the freezeout temperature. This 
technique ceases to be valid near resonance and when the Breit-Wigner expression for the 
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Figure 6: Ratio of the results using interpolation in M — ttidm space and the cross section approxima¬ 
tion detailed in the text to using a full interpolation in M — ttidm — #dm • g q space. The cross section 
approximations are conservative in the bright yellow (light) areas, and overestimates the limit in the dark 
blue (dark) areas. The black dots are interpolation knots in M — ttidm space. Note the ratio takes values 
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are conservative there. See the text for further details. 


width breaks down. We therefore restrict our calculations to regions where Y/M < 1 and 
M > mDM- We use the annihilation rate to quarks for our model as calculated in Ref. [74]. 

We indicate on Figures 2 and 3 a contour, within which the LHC constraint on the 
coupling is stronger than the coupling which would give thermal relic DM. As mentioned 
in the text, for regions inside this line, the coupling strength is constrained to be less than 
the coupling which gives the correct relic density. Hence, the annihilation rate is smaller 
than required, and the relic density will naively be too large. For DM to he in this region, 
either the thermal relic scenario must break down, or the DM annihilates via additional 
channels not considered here. 
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